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We review the prospects of detecting supersymmetric dark matter in the frame- 
work of the Constrained Minimal Supersymmetric Standard Model, and compare 
indirect with direct detection capabilities. 



Recently, both theoretical considerations and and a wealth of experi- 
mental data in cosmology have converged towards a AC DM flat and black 
universe, with the following amounts of dark energy and cold dark mat- 
ter: f2A ~ 0.7, ^CDM ~ 0.3. This last fraction could be incarnated by a 
bath of Weakly Interacting Massive Particles (WIMPs), whose annihilation 
stopped when the universe expansion separated them enough from each 
other, leaving a non relativistic relic density. In the MSSM framework, as- 
suming R-parity conservation, the Lightest Supersymmetric Particle (LSP) 
is stable and is the lightest neutralino (= the neutralino(s) x) in most re- 
gions of the parameter space. If present in galactic halos, relic neutralinos 
must accumulate in astrophysical bodies (of mass Mb) like the Earth or 
most importantly the Sun^, which then play the role of cosmic storage 
rings for neutralinos. The capture rate C depends on the neutralino-quark 
elastic cross section; cr^-g- Neutralinos being Majorana particles, their 
vectorial interaction vanishes and the allowed interactions are scalar (via 
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Xq ^^4^ XQ i'^ ^ channel and XQ A in s channel) and axial (via XQ XI 
in t channel and xg -2^. in s channel). Depending on the spin content of 
the nuclei A'' present in the body, scalar and/or axial interactions are in- 
volved. Roughly, C - ^ X^AT ^bfN J'^ZiM < '"e'sc >^ F{v^,Vesc-, m^, mN), 

where p^, are the local neutralino density and velocity, /jv is the density 
of nucleus N in the body, the nucleus-neutralino elastic cross section, 
Vesc the escape velocity and F a suppression factor depending on masses 
and velocity mismatching. Considering that the population of captured 
neutralinos has a velocity lower than the escape velocity, and therefore 
neglecting evaporation, the total number of neutralinos in a massive 
astrophysical object depends on the balance between capture and annihila- 
tion rates: N^^ = C—CaN^, where Ca is the total annihilation cross section 
(^x-x ti^ss tlis relative velocity divided by the volume. The annihilation 
rate at a given time t is then: 

= IcaN^ = J tanh^t^Ca^) (1) 

with ~ § = csfe when the neutralino population has reached equi- 
librium, and « ^C^CAt^ in the initial collection period (relevant in 
the Earth). So, when accretion is efficient, the annihilation rate does not 
depend on annihilation processes but follows the capture rate C and thus 
the neutralino-quark elastic cross section. The neutrino differential flux 
resulting from xx annihilation is given by: 

Ta fdN\ 

where R is the distance between the source and the detector. Bp is the 
branching ratio of annihilation channel F and {dN/dE)F its differential 
neutrino spectrum. As the direct neutrino production xX ~* exactly 
vanishes in the massless neutrino limit, neutrino fluxes mainly come from 
decays of primary annihilation products, with a mean energy E,^ ^bt 
to ^ (see figure 1). The most energetic "hard" spectra come from neu- 
tralino annihilations into WW or ZZ, and the less energetic "soft" ones 
come from hh. Neutrino telescopes use the Earth as a target for converting 
the muon component of these neutrino fiuxes into measurable muons (see 
S. Cartwright, these proceedings). As both the chargcd-current cross 
section on Earth nuclei and the produced muon range are proportional to 
high energy neutrinos are easier to detect. 

The branching fractions Bp are thus relevant neutralino properties that 
depend on the particular SUSY model considered. We^ have studied these 
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Figure 1. Energy dependence of secondary neutrino fluxes (in arbitrary units) from 
dominant neutralino annihilation channels, as functions of the neutrino energy in units 
of the neutrahno mass. 



in the Constrained Minimal Supersymmetric Standard Model (CMSSM, 
a.k.a. mSugra), whose attractiveness comes from a tractable number of free 
parameters: mo (common scalar mass), mi ^2 (common gaugino mass), Aq 
(common trilinear term) and sign(fi) (supersymmetric scalar mass term), 
all fixed at a high energy scale Equt ~ 2.10^^ GeV, as well as tan/3, fixed 
at the EW scale. The coexistence of these widely different scales introduces 
theoretical uncertainties on the exact definition of the model (especially at 
large tan/3), but the advent of more and more reliable Renormalization 
Group Equations codes (like Suspect2 . 005^ used in this work) tends to 
reduce these. As a bonus, coping with RGE's from the start guarantees the 
expandability of the model to high energies which is the main motivation 
for introducing SUSY and neutralinos in the first place. 

As seen on figure 2, the hard spectra from W^W~ and ti, are found 
at large rng for fixed mi/2 larger than the corresponding threshold. In this 
"focus point" region**, the neutralino has a sizeable higgsino component 
hfracix'^) which allows its annihilation into gauge bosons via i-channel 
gaugino exchange, with a cross-section aA oc ^/rac(x'')''-/rac(x^) ^-^d an 
interesting relic density may survive. Although this region seems a small 
fine-tuned corner of the (mo,TOi/2) plane, relaxing universality may help 
in this respect^'^ Otherwise, the neutralino is an almost pure bino mainly 
annihilating into bb through s-channel A exchange or i-channel sfermion 
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2. The dominant neutralino anniliilation brandling ratios for a typical large 
45 case, as functions of mo (common scalar mass) and mi/2 (common gaugino 



exchange; a low enough relic density can only be found at small mo for 
fixed but not too large mi/2- 

To study the detectability of mSugra dark matter, we have used the 
DarkSusy'' code and computed 1) the relic density, 2) the solar muon flux 
seen by neutrino telescopes and 3) the scalar elastic cross section a^J^p 
relevant to Germanium or Xenon direct detection, for a wide range of such 
mSugra models: mi/2 G (50, 1000) GeV, mo S (0, 3000) GeV, tanP = 
10,20,35,45,50, fj, > 0, Aq = -800,-400,0,400,800 GeV (for tan/3 = 
20,35 only). Among these, we kept only those satisfying the following 
accelerator constraints: BR{b S7) G (2.2, 5.2) x lO"'', a^"^?' G (-6, 58) x 
10-10, m^+ > 104 GeV, ruh > 113 GeV. 

Xi 

In left figure 3, these models are sorted according to their neutralino 
mass and the muon flux $^ above 5GeV originating from the Sun. This 
is compared with past and future experimental sensitivities assuming the 
hardest neutrino spectrum of figure 1 normalized to $^: the lower threshold 
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Figure 3. Indirect (left) and direct (right) detection potential of present and future 
experiments for mSugra models satisfying present experimental limits and offering inter- 
esting relic densities. 

of Baksan or SuperK thus show as a better sensitivity at low m^. Applying 
the shown cuts on the relic density'' separates the models in the 2 rough 
classes indicated above: the lower half corresponds to binos, while the up- 
per half is populated by models in the "focus point" region and neutralinos 
with a sizeable hjradx)- this region, one clearly notices the W^W~ 
and ti thresholds at to^ — 89 and 175 GeV respectively. Between these, 
the neutrino spectrum is indeed hard, and we see that Antares has the po- 
tential of detecting the models with the expected relic density = 0.3. For 
fixed TO^, one also notices the correlation $^ oc which can be un- 

derstood as both the annihilation amplitude (determining the relic density) 
and the spin dependent collision amplitude (determining the capture in the 
Sun and thus the muon flux) are oc h'j^^^x)- When the mSugra neutralino 
is a bino, its spin dependent capture in the Sun is much reduced and the 
muon flux is far below present or future detection abilities. Similarly, neu- 
tralinos captured and annihilating in the Earth give far too low fluxes for 
mSugra models. 

Turning to direct detection, the right figure 3 shows that for small 
masses, both the bino and focus region neutralinos are within reach of 
the next generation of direct detection experiments like Edeweiss II. The 
smaller vertical spread can be traced to the fact that the spin independent 
(or scalar) collision amplitude is proportional to only one power of /i/rac(x): 



The dimensionless Hubble parameter squared is about 0.5 
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Figure 4. Other projections of the same mSugra models as in fig 3: left in the {rrix , ) 
plane; right in the (crj'=_f^"'', plane for Uh^ e (0.03,0.3). Each model is labelled 
according to its detectability: by Antares only, by Edelweiss II only, by both or none. 

which results in a much weaker increase in the focus region. Also notice 
that an ultimate direct detection tool like Zeplin in a maximal version seems 
to cover all interesting relic densities. However coannihilations with staus 
(not included here) can allow for larger masses which would still be out of 
reach. 

Another way to compare the merits of (in)direct detectors of mSugra 
dark matter is shown on figure 4. On the left, all mSugra models of the 
set defined above are placed in the (mx,r2/i^) plane and sorted by their 
detectability. On the right, the models with a relic density fth^ G [0.03, 0.3] 
are placed in the (aj'^p"'', plane and sorted the same. Notice again the 
split in 2 groups, the upper half one again being the mixed neutralinos of 
the focus region. A complementarity between direct and indirect detection 
emerges from this splitting. 

References 

1. G. Jungman, M. Kamionkowski, and K. Griest. Phys. Kept., 267:195-373, 
1996. 

2. V. Bertin, E. Nezri, and J. Orloff. Eur. Phys. J., 026:111-124, 2002. 

3. A. Djouadi, J.L. Kneur, and G. Moultaka. hep-ph/0211331. 

4. J. L. Feng, K. T. Matchev, and F. Wilczek. Phys. Lett, B482:388-399, 2000. 

5. V. Bertin, E. Nezri, and J. Orloff. hep-ph/0210034. 

6. J. R. Ellis, Keith A. Olive, and Yudi Santoso. Phys. Lett, B539:1Q7-118, 2002. 

7. P. Gondolo, J. Edsjo, L. Bergstrom, P. UUio, and T. Baltz. Darksusy program, 
http://www.physto.se/ edsjo/darksusy/. 



